Frustration and sound attenuation in structural glasses 
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Three classes of harmonic disorder systems (Lennard- Jones like glasses, percolators above threshold, 
and spring disordered lattices) have been numerically investigated in order to clarify the effect of 
different types of disorder on the mechanism of high frequency sound attenuation. We introduce 
the concept of frustration in structural glasses as a measure of the internal stress, and find a strong 
correlation between the degree of frustration and the exponent a that characterizes the momentum 
dependence of the sound attenuation r(Q)~Q". In particular, a decreases from in low- 

frustration systems (where d is the spectral dimension), to ~2 for high frustration systems like the 
realistic glasses examined. 



PACS numbers; 63.50.+X, 61.43.-j, 61.43.Fs 

The nature of collective excitations in disordered solids 
has been one of the major problems of condensed matter 
physics during the last decades; the recent devolopment 
of a high-resolution inelastic X-ray scattering (IXS) facil- 
ity Q made Brillouin-like experiments possible in the re- 
gion of mesoscopic exchanged momenta Q~1^10 nm~^, 
which led to the realization that propagating sound-like 
excitations exist in glasses up to the Terahertz frequency 
region. The quantity that characterizes the collective ex- 
citations, and which has been determined experimentally 
in many glasses and liquids is the dynamic struc- 
ture factor S{Q,Lu). Although specific quantitative dif- 
ferences exist among systems, the following qualitative 
characteristics are common to all the investigated mate- 
rials: i) there exist propagating acoustic-like excitations 
for Q values up to Qm, with Qm/Qo^O. 1-^0.5 (where Qo 
is the position of the maximum in the static structure fac- 
tor), which show up as more-or-less well defined Brillouin 
peaks at fl{Q) in S{Q,uj), the specific value of Qm/Qo is 
correlated with the fragility of the glass; ii) the slope of 
the (almost) linear Q{Q) vs Q dispersion relation in the 
Q—i-0 limit extrapolates to the macroscopic sound veloc- 
ity; Hi) the width of the Brillouin peaks, r{Q), follows 
a power law, T{Q)=DQ°', with a«2 whithin the statis- 
tical uncertainties; iv) the value of D does not depend 
significantly on temperature, indicating that the broad- 
ening (i.e. the sound attenuation) in the high frequency 
region does not have a dynamic origin, but that it is due 
to disorder Q. These general features of S{Q,uj) have 
been confirmed by numerical calculations on simulated 
glasses obtained within the framework of the mode 
coupling theory and, more recently, ascribed to a 
relaxation process associated to the topological disorder 
|7|. However, except for the simple 1-dimensional case 
|B|, the widespread finding T{Q)=DQ^ , has not yet been 
explained on a microscopic basis. 

To this end in this Letter we investigate, in the har- 
monic approximation, systems showing disorder of differ- 
ent characteristics, and the role played by the latter in 
determining the value of the exponent a. The analyzed 



systems show either substitutional (bond percolators and 
spring disordered systems on a lattice) or topological 
disorder (model glasses obtained by molecular dynam- 
ics (MD) simulations). We also introduce the concept of 
"frustration" both for structural glasses and, under ap- 
proppriate conditions, for lattice-based systems; we give 
a measure of it, and find a correlation between the value 
of a and the degree of frustration. In particular, for all 
kinds of disorder we find that when frustration is absent 
Q!«4, while increasing the frustration lowers the value of 
a towards and even below the value 2. 

Our samples fall into three major classes: bond per- 
colators above threshold with identical springs connect- 
ing nearest-neighbor occupied sites, topologically ordered 
systems with spring disorder, and model glasses obtained 
by MD. The first two classes are based on a cubic lat- 
tice, while the third represents topological disorder. In 
the first class the disorder is due to the presence of un- 
connected or void sites of the lattice Q; in the second 
class, recently studied by Schirmacher et al. all the 
sites are occupied and nearest-neighbor atoms interact 
through springs whose elastic constants are randomly 
chosen from a gaussian distribution which may include 
or not a negative tail. In both classes the elasticity is 
scalar and each atom is allowed only one degree of free- 
dom. The atoms of the topologically disordered class 
interact through variants of the 12-6 Lennard-Jones (LJ) 
pair potential VLj(r)=4e[((T/r)^^ — (tr/r)^], i.e.: (a) stan- 
dard LJ, V^'^\r)—VL.j{r)] (b) LJ with interaction cut at 
the inflection radius. To = (26/7)i/V, v''''\r)=VLj{r) 
for r<ro and V'^^\r)—{) for r>ro, so that practically 
only nearest neighbors interact; (c) hJ modified to al- 
ways have a positive curvature, V^"^^ {r)—VLj for r<ro 
and V^''^r)=2VLj{ro)-VLj{r) for r>ro; (d) LJ without 
the actractive part (soft sphere): V^'^\r)—Ae{a/rY'^; (e) 
binary LJ mixture of 80% atoms of type A, 20% of type 
B, having the same mass but different LJ potential pa- 
rameters for AA, AB, BB interactions The LJ-like 
glasses were obtained by a microconical simulation which 
is performed at high temperature and followed by a slow 
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cooling down to a state close to the melting temperature 
T,„. From this state a fast quench to low temperaure 
(r/r™«0.1) is performed. After a MD run at low T (to 
ensure that the system is far from the crystalline state) 
the minimum (glassy) configuration ({a;°(z)}i=i..Ar) is 
reached by the steepest descent method. 

For each LJ-like potential, two sets of glasses were pre- 
pared: one had iV=2000 atoms and its dynamical ma- 
trix was diagonalized, while the other, 10000 atoms, was 
used for the calculation of S{Q, uj) by the method of mo- 
ments . The bond percolators had dimension 80^ and 
two bond concentrations, 0.45 and 0.65 respectively, the 
spring-disordered lattices had dimension 80"^; for both 
S{Q, Lo) was calculated by the method of moments. In the 
latter systems, the spring constants k were extracted by 
a gaussian distribution with unitary mean and variance; 
the distribution was cut on the low-A: side at fco=0.3, 0, 
-0.6, and -1. Runs on elongated systems (20x20x3000) 
were also performed, reproducing the results of the cu- 
bic samples in the common Q region. In all cases, the 
width of the Brillouin peaks in S{Q,oj) was estimated 
both by a lorentzian fit and by evaluating the width of 
the spectrum at half height; even when -at high Q- the 
S{Q,U!) does not resemble a lorentzian too closely, the 
two procedures gave results in excellent agreement. 

All the LJ-like systems have the characteristic that al- 
though the system is in a minimum of the total potential 
energy and is stable, all pairs of interacting atoms are not 
at equilibrium distance, which produces internal stress in 
the system It is also possible that some of the par- 
ticles are -at equilibrium- in a maximum of the potential 
energy; these situaions also contribute to the presence of 
stresses in the glass Let V{i) be the potential en- 

ergy of the i-th particle at T—0, and let a single particle 
be displaced from its equilibrium position by an exter- 
nal agent; then all surrounding particles are no longer at 
equilibrium and, if the system can relax, they will move 
towards the new equilibrium position, where each particle 
in general has a different value of potential energy, V'{i); 
this can be either smaller or larger than V{i). We define 
as "frustrated" those particles which have V'{i)<V{i), 
i. e. those particles that, under a small external per- 
turbation, can relax towards a more energetically com- 
fortable situation. The "external perturbation" can be a 
normal mode itself; under the effect of a single normal 
mode, all particles are displaced from their equilibrium 
positions along the direction indicated by the eigenvec- 
tor and, therefore, they change their potential energy; as 
we will see, the change in potential energy can be nega- 
tive. More precisely, when the p-th normal mode -with 
eigenfrequency ujp and eigenvector ep{i)- is switched on 
in the system, and therefore the particle positions be- 
come x° (i) + aep{i) , the i-th particle changes its potential 
energy by £p{i), which is easily written in terms of the 
system eigenvectors as: 



£p{i) = -May 4 S„;5S, Dff {ep{a,i) - ep{(3,j)f , (1) 

where D is the dynamical matrix. Obviously £p{i) de- 
pends trivially aon the amplitude of the normal mode, a, 
and on other system-dependent quantities like the parti- 
cle mass, M , and the interaction strenght, represented by 
the largest eigenfrequency ujo- The relevant quantity is 
therefore obtained by defining the adimensional quantity 
£p{i)—£p{i) /Ma'^ul. In the case of an ordered structure, 
the eigenvectors are plane waves and £p{i) is independent 
on the particle i. In the disordered structure, on the con- 
trary, for each mode there is a distribution of £p{i). This 
distribution (an example is reported in the inset of Fig. 2 
for a mode of full LJ with uj/uJo = 0.03), characterized 
by its mean value /ip {—uj'^/2uj'^) and variance Up, can be 
so broad that many atoms really decrease their energy 
under the action of a mode and, therefore, they are frus- 
trated. To be quantitative, we can define the "index of 
frustration" , Fp as the standard deviation of the distri- 
bution of £p{i). In order to characterize the disordered 
structure itself with its degree of frustration, we can use 
the low frequency limit of Fp] indeed a normal mode in- 
duces in the sample a stress that, by itself, produces an 
" uncofortable" situation for the atoms. Only in the long 
wavelenght (low frequency) limit the stress induced by 
the normal modes becomes negligible and the frustra- 
tion truly associated to the disordered structure shows 
up. The evaluate Fp we need the normal modes, which 
means diagonalization of the dynamical matrix, and this 
puts an upper limit to the dimension of the investigated 
systems. On the other hand, S{Q,uj) is calculated by 
the method of moments , which allows larger samples 
to be treated. The normalized widths, r{Q)/ujo, of the 
Brillouin peaks in S{Q,uj) of the LJ-like systems are re- 
ported in Fig. 1 as a function of Q/Qo- Data from 2000- 
and 10000-atom systems are reported together in the fig- 
ure, indicating that there are no noticeable size effects 
on T{Q). From Fig. 1 we note that: (i) in the low-Q 
region the law T{Q)(xQ°' represents very well the data 
with a values as reported in the legend (see inset, where 
the different data are mutually shifted and the full lines 
represent the power laws); (ii) there is a progressive de- 
crease of slope in passing from the modified LJ system 
(q!=3.5) to the soft potential (a=1.5). The slopes seem 
to coalesce into two groups, for systems a,d, and e Q!Ri2, 
while in the case of b and c, q;>3. We are now going to 
correlate this kind of behavior with the presence and the 
degree of frustration in the samples. 

The frustration index, Fp, of the five LJ-like samples 
is reported in Fig. 2 for some modes spanning the whole 
frequency spectrum as a function of the normalized mode 
frequency LOpjujo- If we consider the variation of Fp, 
the systems can again be arranged into two goups: the 
nearest- neighbor (&) and the modified LJ (c) in the first 
group, the other three systems (a,(i, and e) in the sec- 
ond one. The behavior of Fp is qualitatively the same 
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for all systems down to Wp/wo^O.lS but below this value 
the two groups diverge: Fp of systems of the first group, 
which show exponents a;>3, decreases steadily, while for 
the other three systems (which have a«2) it tends to- 
wards a constant value, Fo, for uj~^Q. The presence of 
such non vanishing low-frequency (and low-Q) frustration 
seems therefore to be strongly correlated to the mecha- 
nism that, in these systems, produces a~2. 
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FIG. 1. Normalized width of the inelastic peaks of S{Q,lu), 
T{Q)/u)o, for the LJ-likc systems as indicated in the legend, as 
a function of the normalized exchanged momentum, Q/Qo- The 
inset shows in an expanded scale the low-Q portion; here the data 
were mutually shifted in order to show the best-fit straight lines 
with the slopes a reported in the legend. 

The previous conjecture is supported by the results ob- 
tained on the lattice-based systems with elastic constant 
disorder. In this case, since the equilibrium positions of 
the atoms are fixed at the lattice sites and the springs 
are all at rest, the systems 



tration only if there are some negative elastic constants: 
the elastic disorder on its own does not produce internal 
stress and frustration. Therefore, in the Schirmacher-like 
spring disorder on a lattice, we would expect a=4 if fco>0, 
a<4 for A:o<0, and a to decrease decreased to 

produce higher frustration index. This is precisely what 
is observed in Fig. 3 for the four studied cases A:o=0.3, 
0, -0.6, -1. As mentioned, we did also evaluate T{Q) for 
elongated samples of size 20x20x3000 in order to reach 
smaller Q- values; in no case did we observe a crossover 
to a dependence for ko<Q although it cannot 
be excluded that it might actually be observed on larger 
samples. It is worth noting that our numerical results 
are in quantitative agreement with a recent evaluation of 
a performed on similar spring-disordered lattices in Ref. 
Ip^ , by using a perturbation expansion on the disorder 
degree, that gives a=4 and a— >1 in the absence and in 
the presence of unstable modes, respectively. 
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FIG. 3. Q-dependence of the width of S(Q,u)) of the lat- 
tice-based systems with elastic-constant disorder, for different val- 
ues of the minimum value (fco) of the elastic constants distribution. 
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FIG. 2. Frustration index Fp for the LJ-like systems stud- 
ied, as a function of normalized frequency, lo/ujo- Inset: a typical 
distribution of €p{i) for a mode of full LJ with ui/uio = 0.03. 

are intrinsically non-frustrated. There may exist frus- 



Further interesting indications come from the study of 
percolating networks above threshold. Here, like in the 
previous case with fco>0, there is no frustration; more- 
over the elastic constants are all equal and disorder is 
produced only by the presence of unoccupied sites. The 
width T{Q)/ujo is reported as a function of Q/Qo in Fig. 4 
for two different concentrations of bonds; we see that 
for both concentrations there are two well defined slopes, 
ai~A and a2~2.3, separated by a crossover; the crossover 
value Qc becomes smaller at smaller concentration. This 
behavior is well understood by considering the fractal 
nature of percolators: for Q<Qc the vibrations experi- 
ence an almost homogeneous system with spectral dimen- 
sion d=3, while for Q>Qc the system is fractal and the 
dynamics is described by the spectral dimension ci=4/3 
Ip^ , ^ . Thus the data of Fig. 4 indicate that in the whole 
Q range r{Q)(xQ'^'^'^^\ where d is the appropriate spec- 
tral dimension. This is exactly the same phenomenology 
observed for the elastic-constant disordered systems on 
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lattice without frustration. Indeed, when fco>0, in both 
3 and 2 dimensions (the latter not reported here) we ob- 
tain r{Q)G<:Q'^'^'^^\ irrespective of the specific value of ko- 
In these systems, like in the percolators, increasing the 
disorder without introducing frustration does not change 
the value of a, but merely shifts the T{Q) curve to higher 
values. On the other hand, the behavior of the highly 
frustrated L J-like systems appears to be different because 
if we try to extract spectral dimensions from the expo- 
nents of Fig. 1 we obtain unexpectedly low values resting 
in the range d=0.6-^l. Such small values are difhcult to 
reconcile with those extracted from the w-dependence of 
the density of states, p{llj)(xuj'-'^~^'' [ p^|Jl7[| which, for the 
present binary mixture, full LJ and soft-sphere potential, 
vields dw2.4. 
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FIG. 4. Q-dependence of the width of S{Q,lij) of bond per- 
colators with two different bond concentrations above percolation 
threshold. 

Thus, all the previous results and discussion indicate that 
frustration may alter the Q-dependence of T{Q). For 
glasses, frustration means that, even if the system as a 
whole is in a minimum of the potential energy, single 
atoms may not be in an energetically comfortable sit- 
uation, which implies that small rearrangements of the 
local structure (due for instance to a propagating density 
fluctuation) may induce rather large displacements of the 
atoms in question towards energetically more favourable 
configurations. We find a strong correlation between low 
a-values and high frustration. This correlation suggests 
a possible microscopic explanation for the observed 
beahvior of T{Q): in the presence of frustration, the 
propagation of density waves is accompanied by large 
and pseudo-random (i.e. not depending on the carrier 
sinusoidal displacement wave) displacements of the more 
frustrated atoms. This entails the presence, in the eigen- 
vector, of large spatial Fourier components other than 
the one corresponding to the dominant Q, i.e. broaden- 
ing. Such broadening is more important at low Q because 
eigenvectors with low dominant-Q are less broadened by 



the "trivial" Q'^ effect found in non-frustrated systems. 

In conclusion, in this work we have studied the role 
of different kinds of disorder in determining the expo- 
nent a, in order to get information on its microscopic 
origin. We have investigated the effects of positive vs. 
negative elastic constants and of substitutional vs. topo- 
logical disorder. In this systematic study we have intro- 
duced a frustration index for the structural glasses under 
consideration. The results present in this work suggest 
that T{Q) behaves hke g^-^+i) in non-frustrated disor- 
dered systems -like for example all the models on lattice 
without unstable modes- and that the exponent a de- 
creases in the presence of appreciable frustration. Values 
a«2 are recovered for the realistic LJ variants examined. 
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